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Abstract 

Based on a recently established formalism (U. Ebert, J. Stat. Phys. 
82 (1996) 183: [1]) we analyze the diffusive motion of a long polymer 
in a quenched random medium. The medium is modeled by a frozen 
semidilute polymer system. In the framework of standard renormal- 
ization group (RG) theory we present a systematic perturbative ap- 
proach to handle such a many chain system. In contrast to [1] we here 
deal with long range correlated disorder and find an attractive RG 
fixed point. Unlike in polymer statics the semidilute limit here yields 
new nontrivial power laws for dynamic quantities. The exponents are 
intermediate between the Rouse and reptation results. An explicit one 
loop calculation for the center of mass motion is given. 
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1 Introduction 

The diffusion of a polymer chain in a quenched environment of obstacles still 
poses many open questions. It is of direct relevance for the motion of a chain 
through a gel. Theoretically two opposit approaches have been pursued: The 
reptation picture in its original form treats the environment as an ordered 
lattice of impenetrable obstacles [2,3]. It neglects all interaction effects ex- 
cept for the impenetrability and it ignores the disordered character of the 
environment, which may lead to entropic trapping [4]. In the asymptotic 
limit of long chains it predicts simple power laws. The opposite approach 
models the environment as a Gaussian distributed random potential with 
correlations of microscopic range [1,5-7]. It thus concentrates on energetic 
effects, and, in a coarse grained sense, also on entropic traps. It however 
completely neglects the impenetrability of the obstacles and it ignores all 
correlations in the quenched environment. This approach results in slowing 
down of the chain motion which is much stronger than predicted by reptation 
theory. In particular it has been argued that the diffusion coefficient of the 
chain decreases exponentially with increasing chain length [1,6]. 



The present work builds upon the random potential model, but includes 
realistic correlations of the obstacles. The environment is modeled as a 
quenched equilibrium configuration of a polymer solution, where the chains 
forming the environment may chemically differ from the mobile chain. In 
the semidilute limit of strongly overlapping background chains this can be 
considered as a good model for a gel, except that we ignore the dynamics 
of short strands of the gel. More serious is the fact that we treat all chains 
as phantom chains. They interact by excluded volume interaction, but of 
course the mobile chain can penetrate the strands of the gel. 

Our model introduces the macroscopic correlation length of the den- 
sity correlations of the environment, which can be identified with the blob 
size £ c (concentration blobs) of the quenched semidilute solution of the back- 
ground chains [8]. Two other important length scales are the radius of gy- 

( B) 

ration R g of the mobile chain and the average radius of gyration R g of the 
background chains. This is in contrast to Gaussian disorder models, where 
the correlation length of the disorder is of the order of the segment size of 
the mobile chain, and where no scale corresponding to R g B ^ shows up. 

Our model turns out to be renormalizable, most similar to standard 
ternary polymer solutions. Renormalization drives the system into fixed 
points, and in limiting situations power law scaling results. The exponents 
are intermediate between Rouse type behaviour and reptation results. 

We performed an explicit one loop calculation of the center of mass mo- 
tion, 

R 2 (t) := <(R cm (t)-R cm (0)) 2 > , (1.1) 

where the angular brackets denote the average over the stochastic process of 
the test chain motion and the bar stands for averaging over the quenched 
semidilute background. Switching off the interaction with the random back- 
ground we simply deal with the free draining limit, which for the center of 
mass motion yields a diffusive behaviour identical to that of a pure Rouse 
chain: R 2 (t) ~ t/n for all times. Here n denotes the segment number of the 
moving chain. With excluded volume interaction between the moving chain 
and the background polymers we find three time regimes: 

(i) As long as the polymer has moved much less than the blob size of the 
background, it does not notice the hindering due to the frozen chains. We 
obtain approximately free diffusion: R 2 (t) = 6D t, where D ~ 1/n, and call 
this regime the ultra-short time regime. 

(ii) What follows is a short time regime, where the moving polymer interacts 
with the background chains. This leads to a slowing down. We find that in 
the case of strongly overlapping chains there is a universal time regime where 
pure subdiffusive power law behaviour emerges: R 2 (t) ~ t 78 . 

(iii) For very long times one always finds diffusive behaviour, R 2 (t) = 6Dt, 
but the diffusion coefficient is reduced as compared to Rouse behaviour. In 
the semidilute limit the diffusion constant decreases with chain length resp. 
monomer concentration as D ~ c~' 45 n~ L45 . 
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Since these exponents are obtained from a first order e-expansion setting 
e = 4 — d — 1 in three dimensions, their quantitative precision is limited. 
Due to the technical complexity of the calculations a higher loop evaluation 
however seems out of reach. 

This paper is organized as follows: In sect. 2 we formulate the model 
and set up the dynamic generating functional which in principle allows for 
a perturbative calculation of any correlation function concerning the mobile 
chain. In sect. 3 we perform a bare one loop calculation of the center of 
mass motion. Sect. 4 is devoted to the renormalization of the model. To 
renormalize the static parameters of the model we establish a connection 
to ternary polymer solutions. The dynamic renormalization is examined by 
means of the one loop result for the center of mass motion. The final scaling 
results for the center of mass motion are presented in sect. 5, followed by a 
concluding discussion in sect. 6. Two appendices are devoted to the relation 
between quenched and annealed averages in polymer statics and discuss the 
case of Gaussian polymers, respectively. 



2 Formalism 
2.1 The model 

We employ the standard Langevin-dynamics of a discrete chain consisting of 
n segments: 

l r ' w = 7 (-a4) H " + f - w )- < 2 ^ 

where !■«(£) e IR ' (i = 1, . . . ,n) fixes the position of the ith segment of the 
moving chain in a (/-dimensional space and d/dr stands for the gradient. The 
thermal random forces fifj.it) (// = 1, . . . , d) are modeled by Gaussian white 
noise. Their second moment is connected to the microscopic mobility 7 of a 
segment: 

(Mt)f jv (lf)) = -5^6(1 - t') . (2.2) 

This guarantees the correct equilibrium distribution V eq , ~ e~ WjV/ for the 
chain (we set /c#T = 1). The Hamiltonian 7i M consists of three parts: 

H M {r 1 (t),...,r n (t);{r {B) }] = 
i 5 ^(r i (0-r i _ 1 (t)) 2 + « o Z d £ S^-v^+Vol^P^i^))- 

i=2 l<«<i<« »=1 

(2.3) 

The first two parts describe chain connectivity resp. excluded volume in- 
teraction in a discrete spring and bead model, where the effective segment 
size is determined by the microscopic length /. The last term describes the 
interaction with the background. 
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The quenched random background consists of M spring and bead chains 
with lengths ni, . . . ,um- A configuration of the background is given by all 
segment positions of all chains: 

( r (BU_ r (B) (B) (B) (B). (B) (B) f ^ 

V v ' V v ' * v ' 

chain 1 chain 2 chain M 

The segment density of the background at point r then reads 

M n m 
m=l j=l 

i.e. the interaction between the moving chain and the background chains in 
eq. (2.3) is nothing else than the excluded volume interaction among seg- 
ments of the moving chain and segments of the background chains. The 
corresponding interaction constant is v , while u in eq. (2.3) governs the 
intrachain-interaction among segments of the moving chain. Note that the 
r^j are independent of t. 

( B) 

Finally we have to specify the distribution function V M of the back- 
ground polymers. We here use a canonical ensemble of chains, just as if the 
background would be a polymer solution in equilibrium: 



'M 

where 



^[{r (B) }] = -^e-^V, ( 2 . 6 ) 



j M n m 

Jj2 ~ r m,j-l) 2 + U Q l d ^2 ^ d ( r mi,ji ~ r m 2 ,j 2 ) 



m=l j =2 all pairs 

( m l Jl);( m 2 <32) 



(2.7) 



is the spring and bead Hamiltonian of M chains with an excluded volume 
interaction constant among all segments of the background and 



4? = Jv{rW}e 



_ n (B) 



(2.8) 



denotes the partition function. The integration f V{y ( - b ">} means summation 
over all configurations of the background polymers in a rf-dimensional volume 
Q — > oo, i.e. 

M n ™ Jd r (B) 

^ ,fl, }=nn7&- 

m=l j=l V ; 
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The factor (Aixl 2 ) d l 2 is introduced for convenience. While angular brackets 
denote the average over thermal forces, a bar stands in the following for 
averaging any observable O over the frozen background distribution: 

= J V{r^}V ( M ] [{r^}]0[{r^}] . (2.10) 

For explicit calculations it is more convenient to switch to the equivalent 
grand canonical ensemble. Thus all final results are understood to hold in 
the thermodynamic limit of infinitely many background chains in an infinite 
volume Q — > oo, where the segment concentration 

c= ^ m ^ nm (2.11) 

is kept finite. A general polydisperse grand canonical polymer ensemble is 
characterized by the function c^n*- 5 -*), giving the concentration of chains of 
length rv- B > G {1,2,..., oo}. The total chain concentration 

oo 

c p = C > (B) ) ( 2 - 12 ) 

n( s )=l 

is connected to the segment concentration c via 

c = Nc p , (2.13) 
where N denotes the average chain length of the background chains 

N= n (B)£^l. (2 . 14) 

In conclusion we should point out that the model presented here is known to 
include all relevant terms in the RG sense, i.e. contributions from any short- 
ranged modification of the model like three-body interactions vanish in the 
limit of long chains. As stressed before our model however does not fall into 
the universality class of short range Gaussian correlated disorder since the 
quenched background has a macroscopic correlation length. This modifies of 
the model pursued in [1] on a long range scale. 

2.2 The dynamic generating functional 

A convenient mathematical tool to calculate dynamic correlation functions 
is the path integral formalism [1,9]. We introduce the dynamic functional in 
Ito-discretization: 

Z M [{h},{h}] := ( e J ^ l=1 V ' J (2.15) 

(2.16) 
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The time discretized meaning of the functional integral is J d[r, f] = 
J Y\ k nj[d d r j{Tk)d d r j(rk) I (27r) d ]. All time integrals without limits are under- 
stood to run from — oo to +00. Differentiation with respect to the external 
fields {hi(t),hj(t)} yields correlation functions of {rj(i), r j(t)}. 

As in [1] we formally have to regularize the theory by adding to the 
Hamiltonian eq. (2.3) a term 

Hn = , (2.17) 

where 

1 n 

R cm (t) = ~J2 r ^) ( 2 - 18 ) 
i=l 

is the center of mass of the moving polymer. This term serves to locate 
the mobile chain in a finite volume of order Q = L d . It is needed for purely 
technical reasons in intermediate steps of the calculation. We will let L — > 00 
as soon as possible. By means of the identity 



8 xdu \ f ddp -■ 



Sd ( r i)= I TT^P^ ( 2 - 19 ) 



the dynamic functional can be written as 



Z M [{h},{h}] = I d[r,r] e V J 



where 



S = J dt^T 7 f i - ir i (*i + 7 



i=i 

is the free action and 



2rj — rj_i — rj + i Rem 
7 ^ 



(2.20) 
(2.21) 



ft = Uo l d J JldtJ2 (?i(*)p) e' p(rj(t) ~ r ^ )} (2-22) 

arises from the excluded volume interaction among segments of the moving 
chain. The interaction with the background chains is contained in Sj B ^: 

S^=v Q l d J J d d re-^p (B \i) J 1 dtJ2^k(t)p)^ prk{t) • (2-23) 



We here use the common abbreviation 

d d p 



(2.24) 
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We now turn to the disorder average, 

Z M [{h},{h}] = J V{v^}V^[{v^}]Z M . (2.25) 

( B) 

Noting that the background coordinates enter Zm only via S) , we rewrite 
this average in terms of a cumulant expansion: Introducing 

.7(f) = -v l d I [ldtJ2 (f fc (0p) e*P( r ^) " (2.26) 
J k=i 



we obtain 



e -SP = e / d*ra(i)p^ (f ) = e Ef=i 7$? , (2.27) 



where 



Sg> = y tffx . . . tfrj a(f . . . a(rj) p^)(r 1 ) . . . p^)(r 3 ) ° 
= (-v l d ) j [ 7*i... 7^ E / 

(r^Pi) • • • (.•,(/,)!>,) e*^^) + • • • + 



(2tt) V( Pl + . . . + p,) jf> ( Pl , . . . , Pi _0 . (2.28) 

The static j-point density correlation function of the background poly- 
mers is defined as (note translational invariance): 

/j B) (pi,..., Pi _i) = 

J d d h . . . dtrj-ie-W 1 * 1 i Pi-i f i-ip(B)(? 1 ) . ..p(B)^ j _ 1 )p(B)(0) ° . 

(2.29) 

The superscript C stands for the cumulant. Now 

/f) = 7b)(0) = c (2.30) 

is the segment concentration of the background polymers. A short glance at 
eq. (2.28) reveals that = 0. Clearly a constant background potential 
cannot exert any force on the mobile chain. The next term j = 2 in the sum 
eq. (2.27) involves a single independent momentum integration, just as Si, 
eq. (2.22). Writing out these two contributions explicitly, we finally find the 
averaged dynamic generating functional as: 

= 777 7 JrTi f , ~S + if dt J™ =1 (hi®*® + 
Z M [{h},{h}] = / d[r,r] e V / 
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exp | - J [u l d J 7 dt^(f,(t)p)e^( r ^ t ) " r M + 

+\vii 2d i?\p) j 7 2 d tl dt 2 jr (^(top) (fjMpyptew - 

ij=l 



cxp 



oo q(B) 

J2-f\- < 231) 



3=3 

The terms tS/^i > 3, contribute only for higher orders (two loop) of per- 
turbation theory and will not be explicitly considered here. 

Before turning to an explicit calculation of the center of mass motion 
in the next section we should compare our general result with the one of 
[1]. Ignoring Sj B ^,j > 3, we see that our result for the averaged generating 
functional resembles that in [1], except for the occurrence of the momentum- 
dependent density cumulant I^\p). The form of the latter is characteristic 
for the quenched background distribution. If we would have worked with a 

Q 

Gaussian distribution of unconnected beads, p( B )(r)p( B )(r') ~ 5 d (r — r'), 
we would have found I^\p) ~ const., just as in [1]. Formally we can map 
the theory as presented here to a theory with Gaussian random potentials by 
taking in the unrenormalized result eq. (2.31) the limit of infinite monomer 
density, c — > oo. For I^\p) the limiting behaviour can be deduced from 
eq. (3.11), whereas a discussion of j-point density correlation functions can 
be found in [10]: 

J (b) (di d- A c -°° i 1 / u o B) ^ for 3 = 2 (9 32) 

1 3 [ p 1 ,...,p 3 _ 1 ) - | q for J>2 • (2.32) 

Thus, with Vq/u^ =: v, the generating functional eq. (2.31) becomes identi- 
cal to that of [1]. Perturbation theory then proceeds in the two couplings u 
and v only. This result has an appealing interpretation since the limit c — > oo 
lets one think of a dense polymer background, where the density correlation 
length should be microscopic indeed. However one should keep in mind that 
our model cannot describe dense systems properly, since it disregards three 
and higher body interactions etc. The limit c — > oo thus is purely formal. 



3 Center of Mass Motion 

We aim at a renormalized one loop calculation of the center of mass motion, 

R 2 (t) := <(R cm (t)-R cm (0)) 2 > . (3.1) 

Remember that angular brackets denote the average over thermal forces, a 
bar stands for averaging over the quenched polymer background. In the 
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formalism developed in sect. 2 we have 

R\t) = -A q , Z M [{h cm },{0}] , (3.2) 

lo 

where 

Hr) an = --(S(t-T)-6(-T)) (3.3) 

Th 

is independent of the segment index i. Expanding the generating functional 
eq. (2.31) in the couplings u and t> , we write: 



z M [{K m }, {o}] = g cm (q, t) + g cm (q, Oi + 0(2 loop). (3.4) 

The tree approximation Gcm(q,t)o can be found in [1], sect. 4.1. For system 
volume L d — > oo (note that the system size L is called £ in [1]) we obtain: 

£ cm (g,t) = y^fle-^ + ^^^i^^cmr^r)) 

= e"^ rot , (3.5) 

where 

r = 1 (3.6) 

n 

is introduced. Considering the one loop contribution we note that the cen- 
ter of mass motion of a Rouse chain with excluded volume (= our model 
without disorder, i.e. Vq = 0) is independent of the excluded volume inter- 
action. Therefore all terms in an expansion of Zm [{h cm }, {0}], which are 
proportional to a power of uq only, strictly vanish. We thus find: 

£ cm(g?t)l = ~v 2 l 2d J d[r,r] e-^ + '/^ELlhiWc^W) 

/ if>(p) / ^dUdhjr (f,(t 1 )p)(f J -(t 2 )p)e«P( r *(*i)- r i^)) (3.7) 
j p J i,j=i 

Except for the integration over p, which involves the momentum dependent 

( f>\ 

density correlation function \P)i we now have to carry through the same 
calculational steps as in [1], sect. 4.2. Taking over the interim result before 
performing the momentum integration from there yields: 

Z M [{h cm },{0}} = e-^|l + t; 2 / M roty/f ) (p)(pq) 2 y t rorfr(l-r/t) 

J- e-P 2 Dij{r) + 2pqr r + Q{2 bop) I (3 _ g) 
i,j=i > 

Here the limit L — > oo has been carried through. D^t) denotes a segment- 
segment correlation function of the Rouse model: 

DiAt) : = 1 <(*(*) - r,(0)) 2 > = 1 1 d[r, f] - r,(0)) 2 e"^ . ( 3 .9) 



An extensive discussion of this quantity can be found in [1], appendix A. 

As a next step we take the momentum derivative in eq. (3.2) and introduce 
variables y = p 2 j_,z = p 2 , where the perpendicular or parallel direction are 
defined with respect to q. Thus p 2 = y + z and pq = ±^/zq 2 . After 
substituting z by z = Dij (r) (y + z) , we perform the y-integration and obtain: 

J™ dz z d l 2 e-"l { 2 B \^z/D l3 {r)) + 0(21oop) | . 

(3.10) 

We now insert the tree approximation for the two point correlation func- 
tion of the background which takes a form well known from random phase 
approximation [11]: 



r(fl ),„, _ cND p (Nq 2 l 2 ^ 

u ( B) l d cND p (Nq 2 l 2 ) 



*n«) = : + o(i ^ . (3 .n) 



The Debye function averaged over polydispersity reads 



D p (x) = \(p(x)-l + x) . (3.12) 

The function p(x) introduced here is the Laplace transform of the reduced 
chain length distribution p(y) := Nc p (yN)/c p [12]. For a monodisperse en- 
semble, p(y) = 5(y — 1), we have p(x) = e~ x , which in eq. (3.12) yields 
the standard Debye function. An expansion of p(x) for small x reads 
p(x) = 1 — x + p"x 2 /2 + . . ., such that D p (0) = p" exists. 

In our result eq. (3.10) we now take the continuous chain limit I — > 0, 
n — > oo with S := nl 2 = const. This amounts to replacing the sums over 
segments by integrals: l A Y17j=i = Jo dsidsj. We then can make all quantities 
dimensionless by dividing by the appropriate power of S. Denoting e = 4 — d, 
our final unrenormalized one loop result reads: 

v n e ' 2 f l T dr 

2(47r) d/2 r(l + d/2) Jo ~ 
s dsids, /A,(t^ ^ d/2 



R\t) = 2dToMl- — — / -^(l-r/t) 



oo 

" ; > / dzz d ' 2 e-* 



S 2 \ S • 
\ l<cND r (NPz/D„(r))° +o(2lQ0p) 



1 + u ( B) l d cND p (NPz/Dijir)) 



(3.13) 
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4 Renormalization 



Eq. (3.13) exhibits a general feature of perturbation expansions in polymer 
theory: The expansion in coupling constants in fact is an expansion in (cou- 
pling constant) • (segment number) 6 / 2 . This can also be seen in the dynamic 
generating functional by general dimensional arguments (power counting) as 
explained in [1], sect. 3.4. For long chains bare perturbation theory in di- 
mension d < 4 thus breaks down. Now the solution to this problem is well 
known. One can map the polymers with microscopic segment size / on equiv- 
alent ones with greater effective segments of renormalized size Ir. Since the 
total extension of the polymers as measured e.g. by the radius of gyration 
remains constant, the renormalized segment number must decrease. This 
way we obtain a renormalized expansion parameter (ren. coupling constant) 
• (ren. segment number) 6 / 2 , which stays finite even for long chains. 

The construction of such a renormalization group (RG) mapping amounts 
to extracting the dependence of the theory on a change of the microscopic 
scale I and to absorb it into the RG mapping. In the continuous chain 
formulation / — > this microstructure dependence is hidden in divergences 
(e-poles) which occur, when the space dimension reaches the upper critical 
dimension d — 4. In such a dimensionally regularized theory the RG mapping 
is set up by introducing renormalization factors which absorb the e-poles 
of the bare theory. To fix these renormalization factors and thus the RG 
mapping, one has to calculate appropriate observables up to the desired 
loop order. This task here however is greatly simplified by the observation 
that polymer statics is contained in the dynamic theory in the sense that 
one equally could calculate all static quantities in the framework of the full 
dynamic theory. This means that the renormalization of all parameters that 
occur also in polymer statics, i.e. all except for r , can be found by an 
inspection of the static theory. This is done in the next subsection. 

4.1 Static renormalization: Quenched vs. annealed 

Calculating a static observable Og like the radius of gyration, we do not need 
the full dynamic theory, but we can employ the static equilibrium distribution 
instead: 

q ' llM n Z M [{r( B )}] ■ 1 ' 

Note that V eq . is taken for a single realization of the background distribution 
and thus the {r^} are to be considered as parameters, not as proper degrees 
of freedom. Note further that in this equilibrium model the segment coordi- 
nates {r} of the test chain are independent of t. The partition function of 
the test chain takes the form 

Z M [{r^}] = J V{r} e-^M; { r(B) }] , (4.2) 
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(4.3) 

IV The static average 

W) = J V{v}V eq .[{v};{v^}]Os[{r};{v(B)}] (4.4) 

and exhibits the well known problem of quenched static averages, namely 
the occurrence of the disorder dependent partition function Zm in the de- 
nominator of eq. (4.1), which in general makes static averages over quenched 
disorder (i.e. the bar in eq. (4.4)) difficult to perform. Now the well known 
solution to this problem is the observation, that the partition function is self- 
averaging, as long as the system size Q is infinitely large compared to any 
other relevant length scale [13-15]. This argument also applies here, since 
we take the thermodynamic limit Q — > oo with keeping chain lengths and 
hence the background correlation length as well as the radii of gyration, 
R g resp. R g B ^ fixed and finite. We rephrase the self-averaging argument in 
appendix A. Let us now discuss the consequences: 

By means of the self-averaging property, Zm = Z M , we can write eq. (4.4) 

as: 

W) = J V{r} O s [{r}; {r<*>}] e"?M{r}; {r^}] / Z M [{v^)}\ . (4.5) 

Now we have a separate average in numerator and denominator, i.e. an aver- 
age as in the corresponding annealed system. But such a system is nothing 
else than the well known static ternary system of two polymer species in 
solution [16-18]. Indeed, taking into account eq. (2.10) we can write out the 
average in eq. (4.5) as: 

V{r, r^}O s [{r}; {r^}\ e^fW; {r^}] - n{f[{r^}] 

(4.6) 

As usual the partition function Z denotes the integral in eq. (4.6) with Os = 
1. This is to be compared with static averages in annealed ternary systems. 
The chains in a ternary solution of Mi polymers of sort 1 and M 2 polymers 
of sort 2 interact via three couplings Uq® = UQ^;a,b G {1,2}; representing 
the intrachain (a = b) resp. interchain (a ^ b) repulsion. Together with the 
chain connectedness the Hamiltonian then reads: 

Hi 2 [{r«;r (2) }] = 

2 M a n£> 2 

4 E E E(rS - -it,-.) 2 + \ E -J - * / «AV WW • 

o=l m=l j=2 a ,b=l J 

(4.7) 



where 

The Hamiltonian 7i M is given in eq. (2.3), with r^t) 
of Oq then reads 
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where 

M a nt ] 

P (a) W = EE^ r - r S) ( 4 - 8 ) 

m=l j=l 

is the segment density of the polymer species a = 1, 2. Let the tracer chain in 
our original system be of sort 1, i.e. Mi = 1, = n, the background chains 
be of sort 2. Comparing + "^m\ eq. (4.6), with the ternary Hamiltonian 
as above, we indeed find complete agreement. The couplings can be identified 
as the three ternary couplings: 

u = u { q 1] , u {B) = u ( o 2 \v = u { d 2) . (4.9) 

Having established the connection of the static behaviour of our model to an- 
nealed ternary polymer systems, we now can simply adopt the corresponding 
renormalization scheme. The one loop results we quote here can be found 
in [18] (note that we set the parameter b u introduced there to b u = 1 from 
the outset). The mapping of the bare static parameters (/; n, N; u , Uq B \v ) 
to the corresponding renormalized quantities (Ir; %, Nr; u, u^ b \v) takes the 
form: 



nl 2 = n R /^|l-- + C(21oop)} (4.10) 

Nl 2 = N R l 2 R ll — + 0(21oopU (4.11) 

(47r)- d /%,r e = jl + -^ + 0(21oop)j (4.12) 

(A>K)- d / 2 ul B) l-* = M ( B )/-i|l + ^ + 0(21oop)| (4.13) 

(47r)- d /^ r e = vl-/- |l + U + ^ B) + ^ + 0(21oop)| (4.14) 

The renormalization of I, n, Uq resp. /, N, Uq is just as in the corresponding 
binary systems. Only the ternary parameter vq is renormalized by all three 
couplings. 

Now aiming at a renormalization of the center of mass motion a glance 
at eq. (3.13) immediately reveals that any possible divergence in the one 
loop integral could be subtracted only by a renormalization of Tq, since this 
is the only parameter that occurs already in the tree result. All e-poles of 
eqs. (4.10)-(4.14) contribute, if inserted into eq. (3.13), only in 2 loop. This 
was to be expected, since R 2 (t) is a purely dynamic quantity. The center of 
mass motion in one loop thus fixes the one loop renormalization of r only, 
to which we will turn now. 

4.2 Dynamic renormalization 

We now search for possible e-poles in the one loop contribution to eq. (3.13). 
Divergences are to be expected for large frequencies and momenta (UV- 
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divergences), that is to say for short times and short segment separation. We 
therefore introduce appropriate dimensionless time and distance variables: 
t := Tqt/S and y : = (sj — Sj)/S/V2f together with s := Sj/S. In these 
variables, the segment-segment correlation function D^t), eq. (3.9), can be 
written as: 

D v( T ) = y/2? p(y, y + 2s/V2f, 2/V2?) . (4.15) 



S 

The properties of the function F(y, z, A) are discussed in great detail in ap- 
pendix A of [1]. We here simply cite the following two representations, the 
latter being suitable for an analysis of the behaviour for short times and short 
segment separation: 



1 , A z y 2 + z 2 
F(y,z,X) = - + A --- + - 



A V3 A A 2 



EA 7r 2 A .2 /A 2 / 2nky 2nkz\ ti 

^P 6 cos— ^ + cos ^- (4.16) 

k=l ^ ' 

oo 

= \y\ + Yl tf(y + ^ - \y + uX \ + f( z + ^ - \ z + uX \] ■ 

v=~oo 

(4.17) 

The function f(y) is defined as 



-y 2 



/(„) = £_ + „ erf (y), (4.18) 
V 71 " 

where erf(y) denotes the error function. The one loop result eq. (3.13), 
written completely in dimensionless variables then reads: 



v n 



e/2 rT 



2dT t " 2(327r 2 ) d / 4 r(l + d/2) J 
•1 rO--s)/V2 



pi poo 

/ drf- d ^(l-r/T) / dzz d '\ 
Jo Jo 



-z 



! ds f dyF(y,y + 2s/V2J,2/V2J)- 1 - d / 2 H(f,z,s,y) 

Jo J-s/Vz? 

+C(21oop) . (4.19) 
The dimensionless time variable T is 

T -%- < 4 - 20 » 

The function H has its origin in the tree approximation of the density cor- 
relation function: 

Hif , i, s, y) = V ^ N °' ( W^,/^)) (4 21) 

1 4- 7/ >l d rN D (— - I 



14 



In the following we will drop the hats in the integration variables f , z, s and 

y- 

Let us first discuss what happens in the case of [1]. To this end we 
as before take the formal limit c — > oo of infinite segment concentration, 
such that H — > v /uq B \ The divergent part in d — 4 then comes from the 
integration: 

pT pi p(l-s)/V¥F 

/ dTT' d/i ds _ dyF(y,y + 2s/V2^,2/V2^)- 1 - d / 2 (4.22) 

JO Jo J-s/\/2t 

For small r we can replace the right hand integrations by 



JO J-c 



dyFfaoo,™)- 1 -*' 2 , (4.23) 

where F(y, oo, oo) = f(y), which can easily be derived from eq. (4.17), taking 

into account the asymptotic behaviour f(x) x ~£°° \x\. In four dimensions 
eq. (4.23) thus yields two times the finite constant 

POO 

/:= / dyf( y y 3 = 3.587... . (4.24) 
Jo 



The remaining r-integral diverges for d — 4 at the lower bound: 

d T T- d ' A = 1(1 + (9(e)) . (4.25) 



>o e 

All together this ends up in a one loop divergence ~ I/e, which has to be 
absorbed in a renormalization of IV 

Now in our system, where c is small and finite, the long ranged density 
correlation function renders the one loop term finite in four dimensions. This 
is not hard to prove if one takes into account that the Debye function obeys 
(cf. eq. (3.12)) 

D p (x) ™ 2/x , (4.26) 

where the right hand side actually is an upper bound. This means that for 
H, eq. (4.21), the following inequality holds: 

< H(t, z, s, y) < —5 y/5?F(y, y + 2s/ V^, 2/y/*?) . (4.27) 

z 

Hence an additional factor y/r occurs, which alters the power of r in eq. (4.19) 
from r~ d l^ to t 1 / 2- ^ 4 and thus improves the convergence for small values of r. 
Together with F(y, z, A) > F(y, oo, oo) it is easy to show that the right hand 
side of eq. (4.19) exists for 2 < d < 6. To see the role of a finite correlation 
length in this connection more clearly we note that the unrenormalized 
tree approximation of the latter in the semidilute limit reads (cf. eq. (3.11)): 

£ {B)2 = -W\ In/f^H^I-. (4.28) 
1 i° u\ >l d cn 
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Concentrating on the dangerous short time part of the time integral in 
eq. (4.19) we according to eqs. (4. 21), (4. 26) write: 



I 



dTT - d/4 2v l d cnV2¥f(y)/z 

l + 2uf ) l d cnV^f(y)/z 



Vq_ 

(B) 



U 



2V2f(y) V R. 



(B) 



(4.29) 



In case of a microscopic correlation length — > the time integral in d = 4 
would not exist. 

The essence of this analysis is that up to the order of one loop T needs 
no renormalization: 

r = r fi {i + o(2ioo P )} . (4.30) 



4.3 Renormalized center of mass motion 

Now we are prepared to renormalize the one loop result eq. (4.19). Since 
in this loop order no e-poles show up, this procedure amounts to a simple 
replacement of the bare parameters by their renormalized counterparts ac- 
cording to eqs. (4. 10)- (4. 14), (4.30). Note that c p = c/N is not renormalized, 
since neither the volume fl nor the number of polymer chains M in the back- 
ground is altered by renormalization. The reduced chain length distribution 
function (see eq. (3.12)) is a normalized function of the RG invariant ra- 
tio n^/N and thus is invariant under renormalization itself. Denoting the 
renormalized time as 

n R l R 



T R = ^, (4.31) 

R 2 {t) vn 



we obtain: 



e/2 r T R r oo 

1 - — „,'.:,? / dr r-W (1 - t/T r ) J dz z d ' 2 e-* 



2dn R l 2 R T R 2 2 + d / 4 r(l + d/2) 

"(i-s)/v^F 

dy F(y, y + 2s/ v 7 ^, 2/v / 27)" 1 " <i/2 

-s/V2t 



fds f 

JO J-i 



H R (T,z,s,y) + 0(2\oop). (4.32) 
The renormalized function H R reads 

H K (r, z, ., y) = (4 33) 
l + u< B >c R N R D p (!^- 



where we defined a dimensionless renormalized concentration as 

c R = { -^f^l d R c P N R . (4.34) 
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To make use of a result expressed in renormalized parameters we have to 
find the explicit RG-mapping, which connects them with the physical pa- 
rameters of the bare theory. This is a standard procedure of renormalized 
perturbation theory, which is done by an integration of flow equations. The 
latter are obtained by taking the partial derivative —Ind/dln of eqs. (4.10)- 
(4.14), keeping all bare parameters fixed. Let us briefly rephrase some well 
known results concerning the flow of the couplings in ternary polymer sys- 
tems [16-18]: The RG flow of the three couplings has 8 fixed points. In the 
limit of long chains, Ir — > oo, there is only one globally stable (IR-stable) 
fixed point, attracting the flow from all directions. This is the symmetric 
excluded volume fixed point, where all three couplings attain the same fixed 
point value: 

u* = u™* = v* = -e + —e 2 + G(e 3 ) . (4.35) 
4 128 

At this fixed point the RG mapping takes the simple form [18]: 




(4.36) 



where Bi, B 2 denote nonuniversal constants, which depend on the chemical 
microstructure of the two polymer species. The critical exponent v governs 
the size of an isolated excluded volume chain without disorder as measured 
e.g. by the radius of gyration, R g ~ n u : 

All the other fixed points are repulsive in at least one direction. The 
symmetric Gaussian fixed point u = — v — 0, in particular, is globally 
repulsive. To reach such fixed points it needs a careful adjustment of the 
nonuniversal parameters of the system. In the sequel we therefore restrict 
ourselves to the symmetric excluded volume fixed point eq. (4.35), which 
governs the generic behaviour in the limit of long chains. A brief discussion 
of other fixed points can be found in appendix B. We note that we could 
evaluate our theory also in all the crossover regime among the different fixed 
points. Due to the lack of data we however do not pursue this program. 

We should stress that the existence of stable fixed points is necessary 
for the occurrence of pure power laws, but by no means necessary for an 
application of RG methods. In the case of Gaussian disorder [1] there is 
in fact no such fixed point in the coupling, which leads to nonlinear scaling 
behaviour. As we will see in the next section, the occurrence of a stable fixed 
point here in the semidilute limit leads to asymptotic dynamic power laws 
with new exponents. 

Having fixed the RG mapping one final step is to be done, namely to 
choose the renormalized length scale. The expansion parameter v n e ^ 2 r -^> 
vn^ in our renormalized result eq. (4.32) suggests to rely on the standard 
choice tir ~ 0(\), which guarantees a finite expansion parameter. This 
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choice implicitly fixes the renormalized length scale Ir to be of the order of 
the extension of the moving coil, since in tree approximation the radius of gy- 
ration behaves as R g ~ nl 2 — > urI r . We should mention that although this 
choice is always possible, in many cases it is not the best one. It, for instance, 
in the semidilute limit leads to logarithmic singularities, which have to be 
cured by reexponentiating the e-expansion. In favorite cases such problems 
can be avoided by a more sophisticated choice of Ir. Quite generally one 
can say that the best and most physical choice of the renormalized length 
scale is to choose it of the order of the smallest relevant length scale in the 
problem. For the static problem of a single polymer in solution there is only 
one macroscopic length scale R g , thus Ir ~ R g . In case of semidilute polymer 
solutions however the blob size £ c (concentration blobs) is the appropriate 
smallest macroscopic length scale [19], while in the context of polymer dy- 
namics time blobs £ t have been invented successfully [7]. Trying to generalize 
these ideas to our system, where both £ c and £ t show up, it turns out that 
they lead to divergences in the one loop result which can not be handled. 
This could have been expected. The reason is that any blob idea is based on 
the existence of two separated regions: Inside the blob, where interactions 
are present and outside the blob, where we essentially have a "free" theory 
without nontrivial interactions or correlations. It is this separation that does 
not work in our case: Looking at short times, such that the time blob is the 
smallest length scale, £ t < £ c , we have inside of the blob dynamic correlations, 
but on scales larger than £ t the theory is not trivial since the segments still 
can be inside a concentration blob, yielding nontrivial correlations. In the 
Gaussian case [1,7] this problem does not arise since the correlation length of 
the disorder is microscopic: £ c ~ / — > 0. For longer times such that £ c < £ t we 
run into the same problems: Although the interaction is screened on scales 
> £ c , we still have nontrivial dynamic correlations. Aiming at a calculation 
of dynamic quantities these correlations ruin the simple picture of a chain 
being composed of uncorrected concentration blobs as in the static theory 
of semidilute polymer systems. 

We therefore rely on the e-expansion together with the standard choice 
of the renormalized length scale by fixing ur — 1. Furthermore we set the 
couplings to their fixed point value u*. Eq. (4.36) then fixes all renormalized 
quantities that show up on the right hand side of eq. (4.32) as a function 
of the bare parameters n, N, c, and t. The bare parameters themselves then 
occur in certain dimensionless scaling combinations only. The corresponding 
scaling variables Y, S, and T are defined as follows: 

. e , Nn = N(B,Y<" =: y (438) 
n R n \B 1 J 

u*c R = ^—L—u*B d 2 c p N dv N R l - du 

n *= =1 cn d »- A ^^u*(^y U B d =:S (4.39) 
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T* f = ^==f- <«<» 

Regarding Tr one should remember that up to the order of one loop T is 
not renormalized, i.e. T R = T = 7/71, cf. eqs. (3. 6), (4. 30). 



5 Results 

Having finished the previous, more technical section, we here first collect and 
interpret our general results for the center of mass motion: The quantity we 
want to investigate is the normalized center of mass motion 



R 2 (t) 

2dn R l 2 R 



It describes the mean squared center of mass displacement R 2 {t) of a chain 
of length n, divided by the radius of gyration R 2 ~ n R l R ~ B\n 2u of the 
tracer chain without disorder. According to the discussion at the end of the 
last section this function depends on three dimensionless scaling variables as 
defined in eqs. (4.38)-(4.40): 

F = F(S,Y,T). (5.2) 

As a function of these variables the normalized center of mass motion T is 
a universal function, i.e. independent of the precise microstructure of the 
polymers. 

The scaling variables have a simple interpretation. T measures time t on 
the scale of the longest relaxation time ~ n l+2u of a free draining chain in 
the excluded volume limit: 

f ~ t/n 1+2v . (5.3) 

The variable Y is proportional to the ratio of the mean chain length iV of 
the background polymers and the chain length n of the moving polymer: 

Y ~ N/n . (5.4) 

The overlap variable S requires a more detailed discussion: In the case of 
a binary semidilute polymer system, chain length and concentration occur 
only in the dimensionless combination ~ c p Rg B ^ d ~ cN du ~ l [8]. The 
overlap describes the number of chains which are present in the volume 
occupied by one chain. In the dilute limit we have — > 0, while the 
semidilute limit is given by -> 00. Now in our case we are interested in 
the properties of a tracer chain of different species and length, which moves 
in such a background. As observed in [18,20], the appropriate overlap type 
variable at the symmetric excluded volume fixed point in such a ternary 
problem is S ~ c p RgN/n (which is called in [18]). If we divide the 
background chains into blobs of length n (assume N > n, see below), then 
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CpN/n measures the number concentration of these blobs in the background. 
S is thus the number of such blobs which overlap in a volume R d ~ n du 
occupied by the moving chain. This variable ensures that it is irrelevant for 
the moving polymer, wether the background blobs belong to different chains 
or not, as it should be, since the moving chain at a given time does not know 
what happens on scales larger than its own extension. Note that S depends 
on c and n only: 

S ~ cn dv ~ l . (5.5) 

The proportionality constants in eqs. (5.3)-(5.5) depend on the chemical mi- 
crostructure of the polymers (cf. eqs. (4.38)- (4.40)). In the framework of an 
RG analysis they can not be calculated reliably but they have to be taken as 
fit parameters instead. 

Quite generally we expect the center of mass motion to become indepen- 
dent of N if the background chains are much longer than the extension of 
the moving chain, Y 1. This conjecture becomes explicit in the one loop 
function H R , eq. (4.33): 

SYDjY 



H R (S,Y;T,z,s,y) = — j-_ r- , (5.6) 



QYD \ Y 5 l ' 

01 ^Py 1 V2TF(y, y +2s/V*f,2/V2T) J 



where the star signals, that we have set the couplings to their fixed point value 
u*. Note further that we have made the dependence on the scaling variables 
S, Y explicit now. Due to the asymptotic behaviour of the Debyefunction 
D p (x) ~ 2/x, cf. eq. (4.26), we in the limit Y — > oo obtain: 



H R (S, oo; r, z, s, y) = ^M±^W , (5 . 7) 

such that the dependence on Y and thus on N in one loop drops out in- 
deed. More problematic is the opposite limit Y — > 0, since in eq. (5.6) 
this limit would compete with the semidilute limit S — > oo. Now hav- 
ing in mind the discussion on the renormalized length scale at the end of 
the previous section, one may rightly ask wether the choice ur — 1 for 
Y — > is appropriate any longer. Since in the limit N n the length 
scale set by the moving chain is much greater than those of the back- 
ground chains, one would argue for a choice TVr = 1 instead. In eq. (5.6) 
this would lead to a replacement of SY by the overlap of the background 
§( B ) ~ c N du ~ l . The limit Y -> then appears to be unproblematic, result- 
ing in S^D p (Yz/V2^F{y,y + 2s/V2T,2/V2T)) Y =° S^p", cf. eq. (3.12). 
But this is not the whole story, since H R then would be independent of r, 
such that the integral J Q Ti? drT^ d ^ in eq. (4.32) diverges in d — 4 at the lower 
bound (a feature which also arises in the semidilute limit, see below). These 
problems occurring in a RG treatment of the limit of a very long chain in 
a system of moderately long chains are known to show up even in polymer 
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statics [8,18,21]. On a heuristic level it is clear that a semidilute system of 
very short background chains will eventually lead to the universality class 
of short range Gaussian random potentials. As a consequence of the new 
singularities for Y — > 0, a precise analysis of this limit would however require 
a reconsideration of the renormalization of the couplings as well. We will not 
pursue this here and restrict ourselves in the following to n<N, excluding the 
limit Y — > (see also appendix B). 

Calculating the explicit scaling function we proceed with an e-expansion. 
In this context two things must be remembered: Firstly it is important to 
note that the renormalized function H R itself, eq. (4.33), must not be ex- 
panded, since such an expansion would destroy the physical idea of screening. 
The structure of Hr comes from a resummation in the treatment of many 
chain systems and thus has nothing to do with the critical behaviour treated 
by renormalization and e-expansion. Secondly, to extract anomalous power 
laws from an e-expansion we as usual have to exponentiate the perturbative 
result. Inserting eq. (4.35) into eq. (4.32), the fixed point behaviour of the 
center of mass motion in first order e-expansion hence follows as: 



HS,Y,f) 



= exp(-^-/ drr- 1 (1 -r/T) I dzz 2 e - z [ 
I 64 y J J 



ds 



dy F{y, y + 2s/ V^, 2/ v / 27)~ 3 

s/V2t 



H R (S,Y;r,z,s,y) + 0(e 2 )j, (5.8) 
where H* R is given in eq. (5.6). 

5.1 Diffusion constant 

We start with analyzing the diffusion constant D. In a way this is the simplest 
quantity, since here the dependence on time drops out due to the limit t — > oo: 

D := lim (5.9) 

t^oo 2dt 

We normalize this quantity by division by the diffusion constant of the free 
Rouse chain, D = 7/71 = r , eq. (3.6). From eq. (5.8) we obtain 

V(S,Y) := 2- = lim gg^lD = exp ( - ± f°° drr' 1 f°° dz z*e~* 
fdsf 1 V 2T dyF(y,y + 2s/V2^,2/V2^)- 3 H R (S,Y;r,z,s,y) + 0(e 2 )). 

JO J-s/V2t ) 

(5.10) 



In the dilute limit, H R ~ S — > 0, cf. eq. (5.6), we of course recover Rouse 
behaviour, 

lim V(S, Y) = 1. (5.11) 

5^0 
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Recall that the intrachain excluded volume interaction u, which is still present 
also for S — > 0, does not affect the center of mass motion as long as the dis- 
order vanishes. For small S eq. (5.10) can be evaluated as a virial expansion 
in powers of S. 

Much more interesting is the limit of strongly overlapping chains, S 1. 
Since lim^^ H* R = 1 (cf. eq. (5.6)), a glance at eq. (5.10) immediately reveals 
that this limit if taken naively would produce a divergence due to the integral 
J drr^ 1 . This was to be expected from the discussion of sect. 4.2, since here 
we obtain a structure similar to the limit c — > oo in the unrenormalized result, 
which is singular in d — 4. The resulting e-pole in [1] induced a nontrivial 
renormalization of r . Here it translates into a logarithmic divergence in 
S. To extract this divergence we exploit the asymptotic behaviour of the 
integrand in eq. (5.10) for small r (cf. eqs. (4. 23), (4. 26)): 

H* R (S, Y; r, z, s, y) 2 1^01Il . (5 .12) 

With this asymptotic form we can perform the dangerous part of the r in- 
tegral near the lower bound in eq. (5.10) and find the desired logarithmic 
divergence: 

fdrr^ r d z*>e-> f ds f°° dyf{y)^ 2S 'f^ 7 * S ~~8JlnS, 
Jo Jo Jo J-oo l + 2SV2Tf(y)/z 

(5.13) 

where I = 3.587. . . is given in eq. (4.24). Isolating this divergent term in 
eq. (5.10) we obtain 

V(S, Y) = S- tI/8+ ° (e2) TZ(S, Y) , (5.14) 
where the residual term TZ stays finite in the semidilute limit: 

Jim H(S, Y) = exp(.56e + C(e 2 )) = 3 1.75 . (5.15) 

5— »oo 

Note that the right hand side of eq. (5.15) is independent of Y and of poly- 
dispersity. This was to be expected, since in the limit of strongly overlapping 
background chains, the moving polymer effectively feels uncorrelated concen- 
tration blobs of the background only. The chain length of the background 
polymers completely drops out. The result eq. (5.14) shows that the diffusion 
coefficient in the semidilute limit obeys a power law. With eqs. (4. 37), (5. 5) 
and d = 4 — e we find: 

D = D V(S, Y) ^ c -^/S+0(e 2 ) n -l-eI/8 + 0(el) ( g Jg) 

In leading nontrivial order of e-expansion this yields in three dimensions 
(e = 1) 

D ~ c- 45 n- L45 . (5.17) 
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Compared to Rouse behaviour, D ~ n" 1 , we obtain a slowing down, which 
however is not as strong as in the reptation picture, D rept _ ~ rT 2 [2]. One 
concludes that the hindering due the disordered excluded volume interactions 
of a frozen semidilute background is not enough to obtain a scaling behaviour 
as if true topological constraints were present. Also the dependence on the 
monomer concentration is comparatively weak. Reptation, combined with 
simple scaling considerations would yield a much stronger slowing down with 
increasing concentration: D rept . ~ c~ x ' zx n~ 2 . The power -1.31 follows from 
the requirements that both D rept ./D must be a function of S only and 
D rep t. ~ n~ 2 . In the reptation case such a scaling argument however is 
dubious, since it neglects the fact that topological constraints generally define 
a new length scale (entanglement length), which is different from £ c , the size 
of concentration blobs [22]. We close this short discussion of power law 
behaviour with the remark that the numerical values for the exponents we 
have obtained should not be taken too serious, since they are based on a first 
order e-expansion. 

In fig. 1 we have plotted the full dependence of the normalized diffusion 
coefficient on S according to eq. (5.10). We have chosen a monodisperse 
ensemble, i.e. all chains of the background have the same length N. This 
simply amounts to a replacement of D p , eq. (3.12), by the standard Debye 
function D(x) = 2jx 2 {er x — 1+ x). Regarding the chain lengths we consider 
the two cases Y — 1, which holds if n = N and the moving chain is chemi- 
cally identical to the background polymers (Bi = B 2 ), resp. Y = oo, which 
describes the case of a long tracer chain that moves in a frozen background of 
much longer chains, N ^ n. Note that the curves are universal, i.e. indepen- 
dent of microstructure effects. In the semidilute limit the diffusion constant 
becomes independent of Y: Both curves obviously coincide for S ^> 1 and 
follow the asymptotic form 1.75 S~ A5 of eqs. (5. 14), (5. 15). Quite generally 
the dependence on Y is very small, a property also known from static quan- 
tities [18]. Note further that V(S, oo) < V(S, 1). The reason is that keeping 
c, n and thus S fixed, an increasing N means that the chain concentration of 
the background c p = c/N decreases. With rising Y the constant number of 
background segments hence belong to fewer and fewer chains. But the more 
the background segments are correlated due to chain connectivity, the more 
effective is the hindering: The diffusion constant drops with increasing Y . 

5.2 Full time dependence 

As explained in the introduction we in a semidilute background expect three 
time regimes: 

(i) For very short times T<S~ 2 the moving polymer does not notice its en- 
vironment, i.e. we expect essentially free diffusion. In this ultra-short time 
regime the polymer moves a distance much less than a correlation length of 
the background: R(t) <C In [1] this regime does not exist, since the 

correlation length is microscopic: — > 0. 
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(ii) What follows is a short time regime T<1 but TS 2 3> 1, where the polymer 
interacts with the background and slows down. The motion then is subdif- 
fusive. 

(iii) In the long time regime T >> 1, i.e. R(t) ^> R g , we obtain a diffusive 
motion of the polymer coil as a whole, but with reduced diffusion coefficient. 

This supposed general behaviour can be confirmed explicitly in our one 
loop result eq. (5.8). Note that in this section we always assume strongly 
overlapping chains: S 1. 

(i) Ultra-short time behaviour: 
For vanishing f we have 

*n = ! , (5,8, 

T^O 1 

which means that there must be an ultra-short time regime where the cor- 
rection to the leading behaviour eq. (5.18) remains small. Indeed, for times 
T<S~ 2 , even for very large S no divergences in the time integral show up. 
This can be seen from eq. (5.12), which reveals that as long as yfrS is small, 
the one loop correction in eq. (5.8) is small. 

(ii) Short time behaviour 

Leaving the ultra-short time regime we analogous to eq. (5.13) have to extract 
the asymptotic behaviour. Using eq. (5.12), we find 



f rfrr- A- Ids /"W ™f'®/' *V 8U»VTS, 
Jo Jo Jo J-oc 1 + 2SV2Tf(y)/z 

(5.19) 

so that 

1 ^ P } = (VTS) K(S, Y, T) . (5.20) 

In the semidilute limit 7Z(S, Y, T) exists and is independent of N and poly- 
dispersity, as it should be (cf. eq. (5.15)): 

lim TZ(S, Y, f) =: 71(f) . (5.21) 

Note that a time dependence still remains, such that eq. (5.20) even in the 
semidilute limit does not yield a pure power law with respect to time. How- 
ever, taking in eq. (5.21) the limit f — > 0, we find: 

K(0) = exp (.30e + C(e 2 )) = 3 1.35 . (5.22) 

This implies that for large overlap there exists an intermediate regime 

l«r£ 2 ,T<l (5.23) 

i.e. const. c - 2 n 3 +M!-0 <g ^ <g const. n 1+2u , where R 2 (t) approaches a power 
law: 

R 2 (t) ~ c -^/S+0(el) n -l + 0^) t l-eI/16 + 0(el) d=3 ^.45^.78 ^ ^ 
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This time dependence again lies between Rouse behaviour, R 2 (t) ~ t, and 
the short time reptation result, R 2 {t) ~ t' 5 . The trivial dependence on chain 
length ~ n" 1 should survive in all orders of e-expansion, if we can think of 
the chain as being composed of n/n^ time blobs of size £ t , which move inde- 
pendently from each other [7]. The latter condition is fulfilled for 1 <C TS 2 , 
which implies £ c <C £ t , so that the semidilute background can not build up 
correlations between the blobs, 
(iii) Long time behaviour 

To describe the long time behaviour T > 1 we have to do the same con- 
siderations as for the diffusion constant, eqs. (5.13)-(5.15). Eq. (5.20) is not 
appropriate since lim^^ TZ(S, Y, T) does not exist. As above we obtain: 

?^±pi) = S- eI / 8+ °^K(S, Y, f) , (5.25) 

Note that lim f ^ K(S, Y , f ) = TZ(S,Y) of eq. (5.14), such that the long 
time behaviour is purely diffusive: 

R 2 (t) T ^ 2dD tS- tI/8+oie2) TZ(S, Y) = 2dDt . (5.26) 

The diffusion constant D has been discussed in sect. 5.1. 

In fig. 2 we have plotted the universal time dependence of T = 
R 2 (t)/2dn R l 2 R according to eq. (5.8) for three values S = 10,100,1000 of 
the overlap. As before we restrict ourselves to a monodisperse ensemble. 
Since the dependence on Y is expected to be weak (cf. fig. 1), we here show 
the case of identical chains only: Y — 1. All curves start out with a Rouse 
like ultra-short time behaviour, then bend down to anomalous diffusion and 
end up with normal diffusion again, R 2 {t) ~ t, but with reduced diffusion 
coefficient D as compared to the Rouse value D . For large overlap one may 
identify an effective short time power law regime (R 2 (t) ~ t 8 for S = 1000), 
the value of the power depending on the overlap S however. Only in the 
asymptotic limit S — > oo, T — > 0, with TS 2 — > oo, we obtain the true 
universal short time power law of eq. (5.24): R 2 (t) ~ t' 78 . 



6 Conclusions 

We have presented a systematic perturbative approach to polymer dynamics 
in a quenched many chain background on the basis of standard RG methods. 
The model describes the hindering due to the excluded volume repulsion of 
a test chain moving in a frozen semidilute polymer background. Topological 
constraints are not included. We carried through an explicit one loop calcu- 
lation for the center of mass motion and identified the scaling variables. In 
the semidilute limit we have found power law behaviour with new exponents. 
They exhibit a slowing down of the polymer motion with increasing overlap 
which however is not as pronounced as in the reptation picture. 
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Trying to compare our results with experimental measurements we are 
confronted with a very sparse amount of data available. Clearly real exper- 
iments always involve entanglement effects. On the other hand one might 
argue that for low concentrations c (but very long chains n, such that we 
still work with great overlap S ~ cn dl/ ~ 1 ) 1 entanglement might not be that 
important, at least for short times. A further problem is that in real exper- 
iments we have an annealed random background, while our theory is set up 
for a quenched average. Only in case of a gel, where the crosslinks fix the 
network to a large extend, or, at least partially, if we deal with background 
chains that are much longer than the test chain one might advocate for a 
separation of time scales, such that the random background may be taken as 
effectively frozen as compared to the tracer polymer. Nevertheless a quan- 
titative comparison between our theory and real experiments does not seem 
to be very promising. 

We therefore must rely on computer experiments. The problem here is 
that most of the work we are aware of focuses on the simulation of polymer 
melts, aiming at a verification of the reptation predictions. This means fairly 
concentrated systems and in most cases annealed averages, i.e. all the chains 
move simultaneously. Kremer [23] in his early work performed one simulation 
with one mobile test chain in a frozen environment of M = 6 other chains, 
each of length n = 200, but he did not vary chain length and concentration, 
making a quantitative comparison impossible. The one simulation he carried 
through moreover belongs to the concentrated regime, where entanglements 
are not negligible. The same objections hold for many other simulations, a 
review can be found in [22]. 

For a real test of the theory one would need a quenched polymer back- 
ground at very low concentration but with very long chains, such that we 
achieve a high enough overlap S ^> 1 to observe the asymptotic scaling be- 
haviour we predicted. Simulations in this direction would be most welcome. 
For a more comprehensive insight into the dynamics of our model we further- 
more started to calculate the time dependent segment-segment correlation 
function, which exhibits a much richer dynamic behaviour. 
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Appendix 



A Quenched vs. annealed in polymer statics 

As announced in the main text, we here will shortly rephrase the self- 
averaging argument of [14,15]: 

Divide the space into boxes of length a. A state of the chain {r} could 
equally be described by the starting point ri and the internal configuration 
{c}, where Cj = r { — r^i (i = 2, . . . ,n). Splitting the integration over r 1 
in an integration over a box volume and the sum over all boxes, such that 
ri = k + s, k e a ■ Z d , we can write eq. (4.2) as 

Z M {{r<»>}] = J V{c] jf J^e-«« . (A.l) 

Now concentrate on one fixed internal configuration s and {c} of the partition 
function in eq. (A.l) (see fig. 3). Then with respect to disorder average all the 
summands (=boxes in fig. 3) in the infinite (Q = L d — > oo) sum Xlkea Z d are 
statistically independent, if all segments in one box are separated more than 
a background correlation length from all segments of all other boxes. This is 
guaranteed by choosing the box size to be of the order a>0(R g ) + 0(R g B ^), 
taking into account that < R g B \ The law of large numbers then yields 
that with probability one X]kea-z d e~ HM = Xlkea-z d £~ Hm holds. 

To illustrate this argument we in [24] carried out an explicit one loop 
calculation for the two point density correlation function between the tracer 
chain and the background chains. This is a static observable, that depends 
on the interchain coupling already in tree approximation. We on the one 
hand calculated the annealed average via diagrammatic perturbation theory 
of ternary solutions, on the other hand we calculated the quenched average 
via the dynamic functional. Note that in this case we used the dynamic 
functional to calculate a static quantity, since in dynamics the quenched 
average is carried out easily (no problems with the partition function in the 
denominator as in quenched static averages). The lengthy calculations within 
the framework of two completely different calculation schemes yielded two 
results, which for Q = L d — > oo exactly coincide indeed. 

B A brief discussion of other fixed points 

In the main text we exclusively evaluated our theory at the symmetric ex- 
cluded volume fixed point, where all three couplings u,u^ and v coincide 
and are equal to the binary fixed point value u* ^ 0. One may ask what 
comes out if we would have dealt with Gaussian chains, where one or both 
couplings u or vanish. To this end one might set e.g. = in our 
one loop result eq. (4.32). A short glance at eq. (4.33) immediately reveals, 
that in this case the limit of strongly overlapping chains Cr ~ S ^> 1 can 
not be taken, since the one loop correction term simply diverges ~ S. (We 
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always assume v ^ 0, since otherwise the test chain is decoupled from the 
background.) 

This problem is closely related to the breakdown of sreening in ternary 
semidilute solutions. Generally speaking the intrachain coupling u of a test 
polymer in a solution of other chains becomes screened, u — > u, where the 
renormalized screened interaction it reads [18]: 

u(z)=u ^ ? A^^. (B.l) 

1 + uWc r N r D p (J£z) 

Here z denotes a momentum type variable. Being interested in the large scale 
properties we can set z — 0. At the symmetric fixed point we then obtain 
complete screening in the semidilute limit: u = for cr ~ S — > oo. (We 
however note that u — u in the limit Nr — > 0: Short background chains can 
not screen the interactions within a long test chain, a property wich refers to 
the fact that we did not analyze the limit Y ~ N/n — > in the main text.) 
On the other hand, in case of either Gaussian background chains, = 0, 
or a Gaussian test chain, u = 0, or both, we obtain u < for large overlap: 
The test chain collapses. 

The physical picture behind that analysis is that a polymer solution of 
two incompatible species demixes. Incompatibility is signaled by v 2 > uu^ 
and eventually leads to a phase separated equilibrium state: A test chain 
is immediately driven out of the background chains. Equlibrium dynamics 
with a tracer chain moving through a semidilute background is only possible 
for a compatible system. With such a system we are dealing in the main 
text. 

In conclusion we note that the center of mass motion in one loop is inde- 
pendent of u, such that u does not show up in eq. (4.32). For other dynamic 
observables like segment correlations we actually do find a contribution sim- 
ilar to u in a one loop calculation. 
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Figure 1: Normalized diffusion constant V(S, Y) = D/D ~ nD as a function 



of the overlap S 
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plotted for a monodisperse background and Y 



1 resp. Y = oo (dashed). Dotted: The asymptotic power law behaviour 
according to eqs. (5. 14), (5. 15). 
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Figure 2: Full time dependence of the normalized center of mass motion 
T ~ R 2 (t)/n 2l/ as a function of normalized time T ~ t/n 1+2u for three values 
of the overlap (S = 10, 100, 1000) together with the free case (5 = 0). The 
quenched background is monodisperse and Y = 1. 
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Figure 3: Illustration of the self-averaging property of the partition function 
of a single chain in disorder. 
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